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On the Determination of Ejecta Structure and Explosion 4-[M/03 
Asymmetry from the X-ray Knots of Cassiopeia A 

J. Martin Laming^ &: Una Hwang^ 


ABSTRACT 


We present a detailed analysis of Chandra X-ray spectra from individual 
ejecta knots in the supernova remnant Cassiopeia A. The spectra are fitted to 
give the electron temperature Tg and (single) ionization age ngt. These quantities 
are compared with the predictions of self similar hydrodynamic models incorpo- 
rating time dependent ionization and radiation losses, and Coulomb electron-ion 
equilibration behind the reverse shock, for a variety of different ejecta density 
profiles described by a uniform density core and a power law envelope. We find 
that the ejecta close to the “jet” region in the NE, but not actually in the jet 
itself, have a systematically shallower outer envelope than ejecta elsewhere in the 
remnant, and we interpret this as being due to more energy of the initial explo- 
sion being directed in this polar direction as opposed to equatorially. The degree 
of asymmetry we infer is at the low end of that generally modelled in asymmet- 
ric core-collapse simulations, and may be used to rule out highly asymmetric 
explosion models. 


1. Introduction 

It has become reasonably well accepted that the star recorded in the astronomer Flam- 
steed’s 1725 catalogue as 3 Cassiopeiae was indeed the supernova that gave rise to the 
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remnant now known as Cassiopeia A (Ashworth 1980). Flamsteed’s observation in 1680 
makes Cas A one of a handful of historical supernova remnants for which we have a precise 
age. In fact Cas A is by now well studied by astronomers at all wavelengths from the radio 
(it is the brightest source in the northern sky, Reber 1944; Ryle & Smith 1948) to TeV 7-rays 
(Aharonian et al. 2001), allowing determinations with a high degree of confidence of impor- 
tant parameters such as the expansion velocities of various parts of the remnant (Fesen 2001) 
and its distance (Reed et al. 1995). This makes it an attractive “laboratory” for studying 
various physical phenomena, such as the plasma physics connected with particle acceleration 
(Vink & Laming 2003), or nucleosynthesis through observed element abundances (see e.g. 
Vink et al. 2001). 

In this paper we analyze the X-ray emission from small localized knots of material within 
the Cas A ejecta in order to try and constrain the ejecta density profile. The sub-arcsecond 
spatial resolution of the Chandra X-ray Observatory has made it possible to study individual 
structures in the Cas A shell. These structures have angular dimension typically ~ 1", which 
is ~ 5 X 10^® cm at the 3.4 kpc distance of Cas A. Consequently a reverse shock of velocity 
1000 km s“^ (a typical value in our models below) will traverse these structures in around 
20 years, a time significantly shorter than the evolution time of Cas A. Thus they may be fit 
with a single value of the ionization age ngt, and models for the evolution of the SNR may 
be applied to infer the timing of the shock heating in a similar way to the use of nuclear 
reaction rates elsewhere in astrophysics, as a powerful tool for investigating the structure of 
Cas A. We emphasise that our goal is not to try and analyze the complete emission from the 
ejecta of Cas A, but that by isolating the X-ray knots we are focussing on those structures 
in Cas A for which we think we have the best physical understanding. In this paper we treat 
only those knots that appear to be 0-rich, in order to constrain the ejecta density profile and 
possible asymmetries in the initial explosion. A companion paper considers in more detail a 
set of Fe rich knots on the east limb, determining element abundances for these knots, and 
from estimates of the Lagrangian mass coordinate for each knot, makes an initial attempt 
to build up a radial profile of the composition of the supernova ejecta. 


2. The Circumstellar Density Profile 

The forward shock speed in Cas A has been recently measured from the two extant 
Chandra ACIS observations with substantial exposures to be in the range ~ 4000 — 6000 
km s"^ (Delaney & Rudnick 2003), with an average value of 4916 km s"^. This is consistent 
with previous estimates based on the observed expansion of the bright ejecta ring (Vink et 
al. 1998; Koralesky et al. 1998; Gotthelf et al. 2001). The average radius is 2.46 pc, leading 
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to an expansion parameter of 0.654. The fact that the ejecta are still clearly visible in X-rays 
has suggested to various authors that Cas A must be in transition from ejecta dominated to 
Sedov-Taylor behavior, and the expansion parameter of 0.654 appears most consistent with a 
uniform circumstellar medium, for which it should vary between 1 0.4 in these two limits. 

Simple interpretations of the spectra of forward and reverse shocked material also indicate 
that the masses of emitting ejecta and swept up circumstellar gas are similar (Fabian et al. 
1980). 

Based on optical spectra, however, the ejecta of Cas A are believed to be rich in O 
(Chevalier & Kirshner 1979), and if this were the dominant constituent of the ejecta, the mass 
inferred from the thermal bremsstrahlung emission measure would be lower by approximately 
a factor of 4 than if it were of an H or He dominated composition. Vink, Kaastra, & Bleeker 
(1996) were the first authors to make this point with regard to X-ray observations, and 
inferred an ejecta mass of ~ 4M©. Favata et al. (1997) revised this mass estimate to the range 
2 — 4M© with further uncertainty arising from how the spectral component responsible for 
the hard X-ray emission behaves at lower (i.e. thermal) photon energies. Recently Willingale 
et al. (2003) also inferred an ejecta mass of 2.2 M©. These lower ejecta masses imply that the 
Cas A blast wave has swept up a relatively higher mass of circumstellar material, dominated 
in composition by H and He, with N also over abundant relative to solar values. In this case 
we expect Cas A to have evolved to the Sedov-Taylor limit, and the observed expansion factor 
would imply a stellar wind solution for the circumstellar density profile, i.e. p oc r“^. The 
existence of a pre-supernova stellar wind is also suggested by the quasi-stationary flocculi 
(van den Bergh 1971), which are presumed to have formed from stellar wind material. 

Another line of argument pointing to a stellar wind density profile is from the existence 
of ejecta at optical emitting temperatures. The reverse shock encounters the densest ejecta 
at the core-envelope boundary, and these are the ejecta most likely to undergo thermal 
instability during the evolution of a young supernova remnant. In the case of a uniform 
circumstellar medium, for typical SNR parameters, the reverse shock only penetrates the 
ejecta core after a time of order 100 years, by which time the ejecta have expanded to 
such an extent that the radiative cooling time is now much longer than the remnant age. 
However with a stellar wind ambient medium, the initial deceleration of the blast wave sends 
the reverse shock into the ejecta core on a timescale of order a few years, when the ejecta are 
still dense enough to radiate significantly. This point was first made by Chevalier &: Fransson 
(1994), and is illustrated in Tables 1 and 2 where the mass of ejecta (assumed pure oxygen) 
that can cool to 10^ K within 320 years of explosion is given for various ejecta envelope power 
laws. Cas A is a little unusual in that it does indeed have optically emitting ejecta, both 
in knots and in nebulosity, again strongly suggesting an s = 2 or similar ambient density 
profile. The existence of radiatively cooled shocked ejecta may also be relevant to the free-free 
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absorption observed towards the center of Cas A in the radio (Kassim et al. 1995; Anderson 
Sc Rudnick 1995). Explanations in terms of absorption by unshocked ejecta pose difficulties 
in that they require unphysically large masses of ejecta. Radiatively cooled shocked ejecta 
at temperatures of a few thousand K can be much more dense, hence dramatically reducing 
the mass of cool plasma required to give the necessary absorption, while producing negligible 
photoelectric absorption of X-rays. The secular decrease of the radio emission from Cas A 
(Reichart Sc Stephens 2000) has also been interpreted in terms of the adiabatic expansion of 
electrons trapped in a shock wave where the magnetic field is also secularly decreasing, as 
one would expect in a remnant circumstellar wind (Berezhko, Piihlhofer, Sc Volk 2003). 

Weaker arguments in favor of an s = 2 external density profile are discussed in Vink Sc 
Laming (2003). The hard X-ray emission out to ~ 100 keV photon energy has been inter- 
preted as due to non-thermal bremsstrahlung from a population of electrons accelerated by 
lower hybrid waves generated by secondary shocks within the Cas A shell (Laming 2001a,b). 
Bremsstrahlung is a relatively inefficient emission process; much more of the electron energy 
is dissipated as heat by Coulomb collisions between accelerated and ambient electrons. Rela- 
tive to s = 0, an SNR expanding into an s = 2 circumstellar medium expands faster at later 
times, and this extra adiabatic expansion allows more heat to be deposited in the plasma 
without a large increase in temperature. Additionally, the existence of thermal instability in 
certain regions of the ejecta allows a much higher power loss by radiation. Consequently, in 
as far as one is prepared to accept the model of Laming (2001a,b), (and this is discussed at 
some length in Vink Sc Laming 2003), a stellar wind external density profile is much more 
plausible than a uniform ambient density. In the same manner, non thermal bremsstrahlung 
emission models are much easier to construct in more radiative heavy element dominated 
ejecta than in an H or He dominated plasma. More recently Delaney Sc Rudnick (2003) 
encountered difficulties in trying to interpret their measurements of the forward shock veloc- 
ity and radius in terms of models using constant ambient density from Truelove Sc McKee 
(1999). Setting the ratio of the forward and reverse shock radii at approximately 3:2 (Got- 
thelf et al. 2001), the models overpredict the observed free expansion rate (determined from 
outlying optical knots). The inferred ejecta mass is also rather small at IAMq. Moving to 
the s = 2 models described below (see section 3.3), the n = 9 models give the correct ratio 
of the forward and reverse shock ratio, the correct relationship between the forward shock 
and free expansion rates, and a more plausible ejecta mass of 2M©. Looking more closely at 
Figure 2 of Gotthelf et al. (2001), the ratio of forward to reverse shock ratio actually varies, 
being about 1.5 on the west limb, but larger ~ 1.8 — 1.9 on the east limb, suggesting a lower 
value for the ejecta envelope power law in this region. This observation will be interpreted 
in more quantitative detail below. 
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3. The Ejecta Density Profile 
3.1. Are the Ejecta Knots really Knots? 

A significant fraction of the ejecta of Cas A are observed in X-rays to be in knots 
or clumps. Some of these are located very close to the forward shock position, inviting 
the somewhat naive assumption that they are significantly overdense compared to their 
surroundings. We argue below based on a number of lines of reasoning that these knots are 
not significantly overdense, but have high ionization ages due to encountering the reverse 
shock early in the evolution of Cas A. Their apparent positions close to the blast wave may 
possibly be due to efficient particle acceleration and stronger resulting shock compression 
at the forward shock, as modelled by Blondin &: Ellison (2001). However we find no other 
indications in Cas A, either morphological or spectral, of shock compressions greater than 4 
that would then be expected, and we do not pursue this issue further in this paper. 

If these knots did indeed undergo reverse shock passage early in the evolution of Cas A, 
then an important question arises as to how they survived to be observed as knots by Chan- 
dra/ACIS. A number of authors (Wang &: Chevalier 2001; Klein, McKee &: Colella 1994; 
Klein et al. 2003; McKee & Cowie 1975; Poludnenko, Prank, Sc Blackman 2001) have mod- 
elled cloud-shock or knot-shock interactions. Upon entering a higher (lower) density medium, 
the shock decelerates (accelerates). The cloud or knot undergoes a different acceleration upon 
shock passage relative to the ambient plasma, which can give rise to a Kelvin-Helmholtz in- 
stability at the interface. In addition the shock inside the cloud or knot will be split into 
transmitted and reflected shocks each time it encounters the cloud/knot boundary with the 
ambient plasma, giving rise to further Rayleigh-Taylor or Richtmeyer-Meshkov instabilities. 
The net result is the destruction of the cloud or knot on a timescale of a few shock crossing 
times. For knots 1" across in Cas A (5 x 10^® cm) the shock crossing time is 15—30 years 
for a 500—1000 km s“^ reverse shock. The knot destruction time of ~ 50 years is consis- 
tent with observed lifetimes of the optical fast moving knots and quasi-stationary flocculi 
(Thorstensen, Fesen, Sc van den Bergh 2001; Sutherland Sc Dopita 1995; van den Bergh Sc 
Kamper 1985; Kamper Sc van den Bergh 1976). However the X-ray knots of interest here are 
of similar temperature to their surroundings (i.e., they have not cooled to optically emitting 
temperatures), and their apparent survival for ~ 200 — 300 years (for models in Table 1) 
places an upper limit on their density with respect to the surrounding plasma of around a 
factor of 3 (see Table 3 of Klein, McKee Sc Colella 1994). These authors did not consider 
underdense knots, but the simulations of Blondin, Borkowski Sc Reynolds (2001) suggest 
that these would be no less unstable. The X-ray knots do not appear to coincide with the 
fast moving optical knots (FMKs), which do appear to be significantly overdense and are 
currently undergoing shock interaction. We believe the X-ray knots underwent reverse shock 
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passage, and perhaps interacted with secondary shocks, some time early in the evolution of 
Cas A and are now expanding with the rest of the remnant plasma, as in fact appears to be 
the case (Delaney & Rudnick 2003). 

X-ray knots of similar density to the surroundings are also much more consistent with the 
apparent factor of only 3-4 increase in surface brightness compared to the ambient plasma. 
Consider a knot with solar abundance ratios of Si, S, Ar, Ca, etc relative to O compared 
to a pure O knot. At temperatures around 1 keV the Si-Ca is predominantly in He- and 
H-like charge states and emits about 30 times more power per ion (~ 15 times more power 
per unit mass; Summers & McWhirter 1979) than O, which is fully ionized and radiates 
only by thermal bremsstrahlung. A knot composed of 0.75 O and 0.25 Si-Ca by mass will 
thus radiate nearly 5 times more energy than a pure O knot integrated over photon energy. 
When considering just emission in the X-ray bandpass, the knot with heavy elements will 
radiate even more compared with the pure O knot by an extra factor ~ exp {E/k^Te) ci 2, 
where E is the lower limit of the bandpass. This occurs because most of the O thermal 
bremsstrahlung power is radiated outside the X-ray band of interest, but the line emission 
from the heavier elements is not. A knot of dimension ~ 10^^ cm compared with background 
plasma extending over 1 pc in depth needs to be around 30 times more radiative in order 
to be visible as a distinct object, so a density enhancement of around 2, giving 4 times 
more radiation, is all that is required. The knots fitted in Table 1 all have Si/0 abundance 
ratios typically 0.5—1 times the solar value (relative to the most recent value of the solar 
O abundance, Allende Prieto, Lambert, & Asplund 2001; Holweger 2001), and also usually 
have small amounts of Fe which we have not considered here. The knot visibility for minimal 
density enhancement becomes even stronger if the surrounding plasma is in part H/He/N 
dominated composition, as would be expected for circumstellar medium shocked by the blast 
wave or ejecta shocked very early by the reverse shock. 


3.2. Observations and Data Analysis 

For our spectral analysis of ejecta knots in Cas A, we use the 50 ksec observation with the 
Chandra Advanced CCD Imaging Spectrometer from January 2000. The observation uses 
the backside-illuminated CCD S3 in GRADED mode, wherein each CCD event is classified 
by its distribution over the detector pixels before telemetry from the spacecraft. We have 
considered three radial series of knots to the east, northeast at the base of the ejecta jet, 
and north-northwest, as shown in Figure 1. The knots have angular diameters of a few to 
several arcseconds, and are numbered starting closest to the center and progressing radially 
outward, except for the base of the ejecta jet, which is shown in more detail in Figure 2. For 
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each knot, we typically extract a spectrum with several thousand counts prior to background 
subtraction. 

The spectra are accumulated in pulse-invariant channels after correction for the spatial 
variation of the gain across the detector, but because the data were obtained in GRADED 
mode, the full removal of the effects of charge transfer inefficiency was not possible. Ap- 
propriate spatially weighted detector response functions and effective area functions were 
constructed for each spectrum. The background was taken from several regions off the 
source. For knots which generally have high surface brightness relative to the background 
(i.e., most of the knots in this sample), the spectral fits are not highly sensitive to the back- 
ground subtraction. These data were obtained relatively early in the Chandra mission (198 
days after launch) so the effect of the buildup of contaminants on the detector is relatively 
small, but we include a component in the spectral model to account for this excess absorption 
(ACISABS in XSPEC). 

For each knot spectrum, we have fitted simple models for a single temperature plasma 
with a single ionization age. The spectral models used do not include emission from the 
element Ar, so a narrow Gaussian line is included in the model to represent the Ar Hea blend 
near 3.1 keV. In all cases, the extremely strong emission lines indicate that the emission is 
dominated by the ejecta, but an additional continuum component is also required (as first 
noted by Hughes et al. 2000). 

We have considered two sets of assumptions for the source of this continuum: the 
elements H, He, C, and N in their solar abundance ratios, and an 0-rich ejecta plasma in 
which H, He, C, and N, and Ne are absent and O provides the continuum. We favor the 
O continuum fits as they are more straightforwardly compared to the model calculations 
described above and are a more plausible description of the composition of the ejecta (after 
Vink et al. 1996). Both sets of fits give reasonably comparable fits, although in specific 
cases, either might at times give a better value fot the fit statistic than the other. We use 
the O continuum fit results, noting that the values of the fitted temperature and ionization 
age are typically affected by less than 15% and 40% respectively if the solar composition 
continuum is used instead, though the effect is larger in a few cases. A comparison of the 
fits for knots 4 in the E and NNW are shown in Figure 3 and Figure 4, for both 0-rich and 
solar compositions. 

Because of residual uncertainties in the gain of the detector (e.g., from charge transfer 
inefficiency), combined with the known spatially dependent radial velocities of the X-ray 
emitting gas (e.g., Markert et al. 1983), we have allowed the redshift of the model to be 
freely fitted; likewise a Gaussian smoothing function whose width scales with energy was 
also included in the model. 
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Table 1 shows the results of the spectral fits for all three radial series, giving the tem- 
peratures, ionization ages, and abundances of the elements Si and Fe (by number relative 
to O). In most cases, we obtain reasonable values of the fit statistic {x^/dof ~1.5). Poorer 
fits (e.g., in the NNW and for the outer knots in the E) generally arise from the presence 
of Fe K emission in the spectrum that is not well-fitted by the model. The fits are driven 
strongly by the prominent Si line emission, but knots that are very rich in Fe relative to Si 
are described by different parameters than the knots with strong Si emission, particularly in 
having a higher ionization age. For knots that show a mixture of these characteristics, fully 
successful models would have to include components to describe Si and Fe separately. 

In principle, the blast wave also contributes to each of the knot spectra. To verify 
that it is valid to neglect the blast wave contribution in examining these compact regions, 
we have compared the fit result for the well- isolated knot NNWl using the standard off- 
source background compared to a local background region surrounding the knot. The local 
background should incorporate the blast wave contribution so that using this background 
spectrum eflFectively subtracts the blast wave contribution at the knot. We found no signifi- 
cant difference in the temperature (kT = ^eV compared to kT = l.OTlolog keV with 

the off source background), and a poorly constrained ionization age that is consistent with 
the value previously obtained. 


3.3. Ejecta Profile Models 

Our methods for computing the ejecta model follow from calculations in Laming (2001b), 
Laming & Grun (2002) and Laming & Grun (2003), and are summarized in the two appen- 
dices to this paper. We take an analytic approximation to the hydrodynamics for supernova 
ejecta expanding into a remnant stellar wind from Truelove &: McKee (1999), specified more 
thoroughly in Appendix A. Within this framework we calculation the time dependent ioniza- 
tion balance for the elemental composition of the ejecta, the electron and ion temperatures 
allowing only for equilibration between species by Coulomb collisions, and radiative and 
adiabatic expansion power losses, as outlined in Appendix B. 

We have calculated two series of ejecta models, summarized in Tables 2 and 3. Figures 
5 and 6 show the locus of electron temperature, Tg, against ionization age, rigt, for these sets 
of models. We concentrate on these parameters since they are readily determined from fits to 
the data, and the curves of Tg against Ugi are, for reasons discussed further below, surprisingly 
robust against various changes to the model, being mainly dependent on the ejecta density 
profile specified through the envelope power law n. Both sets of models have an ejecta mass 
of Mej = 2Mq. The first has E^i = 3 x 10®^ ergs explosion energy and an ambient density at 
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the current blast wave position in the range 3. 1-3.6 H atoms cm“^, depending on the value 
taken for the blast wave radius (i.e. pr^ = 21 H atoms cm“^ pc^)- These values are chosen 
to match the blast wave velocity and radius measurements of Delaney &: Rudnick (2003) as 
closely as possible for values of n ~ 9 (in the middle of our range). The second set has the 
same ejecta mass but an explosion energy of 2 x 10^^ ergs and circumstellar density reduced 
from before by the same factor 2/3. These models may be scaled keeping £'51 a Mej oc prf 
which leaves the blast wave velocity and radius invariant. The first set of models is chosen to 
match the blast wave velocity and radius found by Delaney & Rudnick (2003), but implies 
a mass of shocked circumstellar material of close to 17 M©, significantly larger than the 
~ 9M© found by Willingale et al. (2003). The emission measure determined from thermal 
bremsstrahlung using the BeppoSAX data studied in Laming (2001a, b) interpreted as coming 
from an H-He dominated plasma is consistent with masses of shocked circumstellar plasma 
in the range 6.7 - 10 M© for densities ahead of the blast at its current location of 3 - 2 H 
atoms cm“® (or equivalent mass) respectively. These densities combined with the blast wave 
radius give masses of 11 - 16.8 M© from simple geometrical considerations for densities of 
2 - 3 H atoms cm“^, and so it is clear that more consistent results derive from the lower 
circumstellar density. Using this however gives a radius for the blast wave that is slightly 
too small. We suspect the Cas A may have undergone an initial period of expansion into a 
tenuous WR wind before hitting the much denser red supergiant wind into which it is now 
expanding. We discuss this further in section 3.4. 

In Tables 2 and 3 we give the predicted blast wave velocity, radius and expansion 
parameter 77, the reverse shock radius and the ratio of the shock radii. This ratio is measured 
to vary between 1.52 and 1.73 with lower values being found in the west limb of Cas A 
(Gotthelf et al. 2001). From the models we can see that this immediately implies n ~ 9 
on the west limb decreasing to 6-8 elsewhere. The following table entries give dynamical 
parameters from the shock solutions, and the final column gives the mass of gas that may 
cool to optical emitting temperatures (~ 10'* K) within 320 years of explosion. This mass 
increases with n, being zero for the n < 7 models. The initial inference of highest n on 
the west limb is appealing since this is where the strongest radio and non-thermal X-ray 
emission is located coincident with the position of the contact discontinuity, (Bleeker et al. 
2001; Gotthelf et al. 2001) implying that this is the location of strongest magnetic field. The 
high mass of radiatively cooled and thus very dense gas in this region will lead to stronger 
Rayleigh-Taylor instability at the contact discontinuity here than elsewhere, and hence a 
higher magnetic field. In what follows we will develop the idea of varying ejecta envelope 
power law with position in the remnant and try to quantify the degree of anisotropy by 
reference to our ejecta models. 

We argued above that only a modest degree of overdensity may exist in the X-ray 
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knots, in order for them to survive against hydrodynamic instabilities following reverse shock 
passage. Here we demonstrate that the clumping of ejecta by a factor up to around 3 does 
not greatly aflPect the dependence of Tg on nj,. Consider ejecta with a reference value of Tigt 
at the current epoch, which encountered the reverse shock at time ^o- If this were overdense 
by a factor a, say, then it would have passed through the reverse shock at a time tc = cxto, 
since the electron density varies as rig oc 1/t^ and riedt oc 1/to for to « t' . The clumped 
ejecta are shocked to a temperature lower by l/o: than the unclumped ejecta, but undergoes 
less cooling by adiabatic expansion. In ballistic expansion temperatures are reduced by 
adiabatic expansion by a factor 

T' {t') ~ T (to) exp ^d?j T (to) exp ~ T (to) . (1) 

Hence clumped ejecta shocked at ato will have a final temperature higher than that for 
unclumped ejecta by a factor a~^ x In Sedov-Taylor expansion r oc t^^^ the 

final exponent 4/3 become 8/9, and clumped ejecta has a final temperature lower than that 
for unclumped ejecta by a factor Additionally, if Ti » Tg so that Tg oc (TiUet)^^^, 

(Laming 2001a), the temperature difference factor varies as -> remembering 

that the fits are sensitive only to T^. Consequently only for large overdensities am will 
significant changes in Tg result. We have verified with numerical calculations that these 
conclusions hold, using formulae given by Sgro (1975) for transmitted and reflected shock 
velocities. 

It has been argued previously (Laming 2001a,b) that the non-thermal hard X-ray con- 
tinuum of Cas A extending out to 100 keV is bremsstrahlung from a population of electrons 
accelerated by secondary shocks propagating within the shell. These shocks arise as the for- 
ward or reverse shocks encounter density contrasts and split into transmitted and reflected 
shocks. How justified are we in assuming that the ejecta knots of interest encounter the re- 
verse shock and then expand self similarly with the rest of the plasma, undergoing no further 
interaction with secondary shocks? Following the discussion and analysis in Vink & Laming 
(2003), we believe that the electron acceleration actually occurs in shocked circumstellar 
plasma in the immediate vicinity of the contact discontinuity, since an instability operating 
in colder shocked ejecta is much less likely to yield the required electron energies. Thus it is 
the blast wave encountering quasi-stationary fiocculi that gives rise to the secondary shocks 
of interest, not the reverse shock. The ejecta knots we study here do in fact have Lagrangian 
mass coordinates that place them well back from the contact discontinuity, i.e. there is a 
lot of mass through which a secondary shock originating at the blast wave would have to 
propagate through to reach them, and we consider the possibility of significant heating by 
secondary shocks unlikely. The reverse shock is currently propagating through even more 
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tenuous plasma than the forward shock, and again is unlikely to produce secondary shocks 
that perturb the ejecta knots. In our view, the ejecta were subject to various hydrodynamic 
instabilities following reverse shock passage and interactions with secondary shocks for a 
time comparable to the time taken for the reverse shock to propagate further into lower 
density ejecta. What we now see as “knots” are those structures that were able to survive 
against these instabilities, due to having little density contrast with their surroundings. 


3.4. Discussion 

Figures 5 and 6 plot the locus of Te against Uet for the two sets of ejecta models, with 
results from those spectral fits for which the degree of freedom is less than 2 to the 

knots superimposed. In both cases the knots from the jet base on the NE limb consistently 
indicate the lowest value for n. In each case NE 3, 4, 6, 7, and 8 clearly indicate n < 6, 
and NE 5 (not plotted) is also consistent with this. NE 9 and 10 indicate steeper envelopes, 
n = 7 — 9 depending on the choice of model and NE 1,2, and 11 suggest n > 12 or n = 6 if 
coming from the envelope portion of the curve. We suspect that this relatively large change 
in n, if real, over small apparent distances is mainly due to projection effects, knots NE 1,2 
and 11 are either in front of or behind NE 3, 4, 6, 7, 8, 9, and 10, which all come from 
the brightest part in this region. Comparing to Delaney & Rudnick (2003), our knots are 
closest to their regions 2-5 from which they measured forward shock velocities in the range 
4361-6520 km s“^ and radii 2.08-2.43 pc. The series of knots on the east limb all indicate 
n ~ 9 or higher, except for knots El and E2 which are the two nearest the projected SNR 
center and are closer to n = 6. Delaney & Rudnick (2003) do not take any measurements 
from this region of the blast wave, the closest being their region 29 some distance to the 
south. The NNW series of knots pose more problems in fitting the Fe K region, but those 
for which adequate fits can be found, NNW 4 indicates n = 6 and NNW 5 and 7 give 
n = 7 — 12, depending on the hydrodynamic model, with n = 7 — 8 being preferred. These 
knots are closest to regions 13 and 14 in Delaney & Rudnick (2003). At this point we give 
more attention to the model with £^51 = 2 and prl = 14. Reasons for doing so are given 
above with reference to the shocked circumstellar material and the geometry of the blast 
wave. We also find that the lowest temperature knots for the higher energy model are only 
consistent with ejecta envelope power laws n = 30 — 50 which are much steeper than any 
realistic explosion model would predict. Hydrogen rich ejecta give electron temperatures 
lower for a given rigt by a factor 2-3. More electrons per baryon means that the equilibrated 
electron temperature will always be a factor of 2 lower than for O. Additionally, the lower 
ion charge leads to slower electron-ion equilibration before this point is reached. We consider 
that these points reinforce our assumption of 0-rich ejecta. 
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We would like to interpret the diflFerent ejecta envelope power laws as due to asymmetries 
in the explosion. The models of Tables 2 and 3 give different ejecta core densities. Varying 
E 51 to keep the core density invariant yields 


E,i oc Mf 


(n - 3)^^^ 

nV3 (ji — ^y 


( 2 ) 


and we can see that the lower values of n actually correspond to higher explosion energy in 
that particular direction. An n = 6 model corresponds to 66% more energy than n = 9. 
The variation still leads to blast wave radii predicted to vary between 2.35 and 2.62 pc, with 
lower n giving larger radii. This is similar to the range of blast wave radii found by Delaney 
& Rudnick (2003), but is not completely consistent with the locations they measure, and so 
for the time being we further adjust all models to give the same blast wave radius. These 
new models are given in Tables 4 and 5. The degree of explosion asymmetry one would infer 
from the shallower ejecta profile is now slightly reduced; n = 6 corresponds to about 40% 
more energy than n = 9 for our preferred hydrodynamic model. The curves of Tg against 
rigt for these new models for n = 5.5 and 6 are plotted as dotted lines on Figures 5 and 6. 
For higher n the new curves are negligibly different to the previous ones. We can see that 
for both hydrodynamic models the NE knots at the jet base now clearly suggest n = 6 and 
not n = 5.5, and other conclusions remain the same. 


The degree of energy asymmetry in the initial explosion is at the lower end of that coming 
from simulations. Two basic mechanisms of asymmetric core collapse explosions have been 
discussed in the literature. Fryer &: Heger (2000) and Fryer & Warren (2003) model the 
explosion of a rotating star in two and three dimensions respectively. Compared to non- 
rotating explosions, the rotating core appears to be stabilized against convection, and the 
core bounce that seeds the neutrino-driven convection is also weakened. Both these effects 
are reduced in the polar regions, and consequently an asymmetric explosion which is stronger 
in the polar regions results. Ejecta velocities about a factor of two higher in polar relative to 
equatorial regions are predicted, corresponding to a factor of four difference in kinetic energy 
if the density is the same in polar and equatorial regions. ®®Fe is generally synthesized (as 
^®Ni) along the axis. Additionally, the generic problem of core-collapse simulation producing 
too much neutron-rich matter appears to be slightly exacerbated, though Fryer Sz Warren 
(2002) caution that much of the input physics upon which this conclusion depends is still 
uncertain. In contrast Khokhlov et al. (1999) model a jet induced explosion. This mechanism 
is independent of details of convection. A magnetorotational instability in the collapsing 
core (recently studied in more detail by Akiyama et al. 2003) accelerates two jets along the 
rotation axis. These drive bowshocks ahead of them that compress most of the remaining 
ejecta into an equatorial torus. Similar anisotropies in ejecta velocities to the rotating case 
above are found, though due to the anisotropic ejecta density distribution the kinetic energies 
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will be less anisotropic. 

Nagataki et al. (1998) postulate similar degrees of asymmetry in their models to explain 
high yields of ^^Ti produced in a-rich freeze out in the polar regions. Their model Al, with 
similar asymmetry to models discussed above, produces 1.8 x lO'^M© of ^^Ti and O.ObM© 
of ^®Ni for a suitable choice of mass cut. The ^^Ti mass is similar to that inferred for Cas 
A (Vink et al. 2001; Vink & Laming 2003). The mass of ^®Ni produced is less reliably 
estimated from observations. In a companion paper (Hwang k, Laming 2003) we do indeed 
find evidence of ^^Ti production in a-rich freeze out, not in the jet region but further out 
along the east limb series of knots. Indeed this, the general distribution of heavy elements 
in a torus around the “jet” axis, and the existence of the “jet” itself in Cas A do resemble 
the model discussed by Khokhlov et al. (1999), though considerably more detailed analysis 
would be required to confirm anything more quantitative. 


The progenitor mass may now be estimated from the parameters of our models. Most 
of the stellar wind is expelled during the red supergiant phase, and it is into this wind that 
we believe Cas A is now expanding. A massive star spends ~ 2 x 10^ years in this stage 
(Garcia-Segura, Langer, k Mac Low 1996), roughly independent of stellar mass. The mass 
loss rate in this phase is given by 


dM 

dt 


= 4irpr^v.u) = 3 X 10 ^ 


/ p^b \ 


V 1 H atom cm'^pc^ / 

VlOO km s"V 


Mq year \ 


( 3 ) 


where is the stellar wind speed and our favored value for prl = 14. We estimate 
from the speeds of the quasi-stationary fiocculi of Cas A. In general, they are observed to 
have proper motions ranging up to around 500 km s~^ (van den Bergh k Kamper 1985), 
consistent with clumps a factor of ~ 10^ — 10^ more dense than the ambient plasma having 
undergone some acceleration during the passage of the blast wave. The lowest observed 
velocity is 20 km s“^ for QSF 10. Taking this value in equation (3) we estimate a total 
mass loss during the red supergiant phase of 17 M©. The iterated blast wave speeds and 
radii in Tables 4 and 5 are respectively slightly higher and lower than the average values in 
Delaney k Rudnick (2003), with our second set of model for lower pr^ being more discrepant. 
This might suggest an early stage of expansion into a much more tenuous stellar wind from a 
Wolf-Rayet progenitor, before the blast wave encounters the much more dense red supergiant 
wind, similar to a model for Cas A proposed by Borkowski et al. (1996). In the case of Cas A 
such a fast tenuous wind can not have existed for very long before explosion because one 
would then lose the radiative instability of shocked ejecta and the accompanying optical 
emission, and so the mass lost during this period is probably negligible. To the estimate 
above we should then add the observed mass of ejecta and compact objects of about 3—4 
Mq and mass loss during the main sequence evolution to get a total progenitor mass in 
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the range 20-25 Mq, at the lower end of Wolf-Rayet progenitor masses inferred by Massey, 
DeGioia- Eastwood, & Waterhouse (2001). The assumed value for the wind speed of 20 km 
s“^ is significantly lower than is usually assumed for red supergiants (e.g. Garcia-Segura, 
Langer, & Mac Low 1996; Lamers & Nugis 2002) where speeds around 100 km s”^ seem 
more plausible. This would increase our mass estimate to close to 100 Mq. This seems to us 
unlikely since the progenitor would then spend much longer in the Wolf-Rayet phase before 
explosion (Woosley, Langer, & Weaver 1993), so much so that the Gas A blast wave should 
still be moving through Wolf-Rayet wind rather than red supergiant wind, and as discussed 
above, no thermal instability should be present in the ejecta. The Wolf-Rayet phase is much 
shorter for less massive stars, again arguing for a progenitor at the lower end of observed 
progenitor masses for Wolf-Rayet stars. These masses are consistent with those estimated 
from the initial ^^Ti mass inferred from 7 ray observations (lyudin et al. 1994, 1997; Vink 
et al. 2001; Vink Sc Laming 2003) compared with explosion calculations (e.g. Timmes et al. 
1996). Rothschild Sc Lingenfelter (2003) provide a more complete discussion of this point. 


4. Conclusions 

This paper has presented an initial analysis of the X-ray ejecta knot spectra from Cas- 
siopeia A. Although considerably more labor intensive and detailed than any previous anal- 
ysis of imaging spectroscopy of supernova remnants known to us, we believe that we are 
close to being able to infer some of the really fundamental aspects of this object. Rather 
than attempting to analyze and model the entire spectrum of the remnant, we have isolated 
spectra from ejecta knots whose properties we think we understand, which when combined 
with 1-D hydrodynamical models allow us to infer ejecta envelope density profile power 
laws and hence explosion energies in different directions in the remnant. The fundamen- 
tal discriminant in this work is the variation of the reverse shock velocity in the models 
with different ejecta density profiles. The faster reverse shock in shallower profiles leads to 
higher ejecta temperatures and a larger separation between the forward and reverse shocks. 
These considerations lead to conclusions concerning the energy asymmetry of the explosion, 
which we find to be around 50% larger in polar regions than at equatorial regions outside 
the so-called jet. This is at the lower range of energy asymmetries coming from existing 
core-collapse explosion simulations, and on this basis alone we are not yet able to distinguish 
between rotating convection driven explosions or the jet induced explosions. Aside from the 
demonstration of new data analysis and interpretation techniques, the one further important 
conclusion we wish to emphasize is that only with significantly deeper Chandra observations 
of Cassiopeia A and other supernova remnants will the full potential of these methods be 
realized. 
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A. Hydrodynamics for s = 2 Supernova Remnants 

Here we summarize the equations governing the evolution of the forward and reverse 
shock velocities and radii, following Truelove & McKee (1999). These authors concentrated 
on the s = 0 uniform density circumstellar medium case, so here we outline the extension 
of their results to s = 2, as implemented by us. The ejecta density profile is taken to be a 
constant density core with an envelope obeying p (r) oc where n > 5. The circumstellar 
medium density profile is p (r) oc where s < 3. For a uniform density ambient medium 
s = 0, and for a steady stellar wind, the case we will take for Cas A, s = 2. We work in 
similar units to Truelove & McKee (1999); 

to = 423M^fF-^^^p-^/^yeeiTS (Al) 

xo = 3.07MVV"^/^PC (A2) 

for s = 0 and 

to= 5633Mff (pR^)"Vears (A3) 

Xo= 40.74Mej (pRir^ pc (A4) 

for s = 2, where Mej is the ejecta mass in solar masses, E 51 is the explosion energy in units 
of 10®^ ergs, p is the circumstellar density at the blast wave in hydrogen atoms (or equivalent 
mass) per cc and Rh is the blast wave radius in pc. The unit of time to is related to the 
so-called Sedov-Taylor time in McKee & Truelove (1995) by tsx = to/2.024. The unit of 
distance is similarly related by a factor 1.377 to their Sedov-Taylor distance xsr- These 
fiducial quantities apply to the case of uniform density ejecta expanding into a uniform 
density ambient medium. For t < tsx the remnant is in the ejecta dominated phase where 
the mass of the ejecta is dominant, and for t > tsx the Sedov-Taylor phase where the swept 
up circumstellar mass dominates the remnant evolution. 

The ejecta dominated phase is further divided into an initial period while the reverse 
shock is propagating through the ejecta envelope, and a later period once it reaches the core. 
During the initial envelope phase of evolution the blast wave radius is given by (this is the 
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generalization of equation (75 of Truelove & McKee 1999, for arbitrary s) 

f 4y> 3(3-s)^ r 10 (n- 5) 

^ ' (!>ed 47t (n — 3) n 3 (n — 3) j 


(A5) 


Here Ied^'^ the ratio of the radii of the forward shock to the reverse shock, known as the “lead 
factor”, so that Rr = Rb/lsD, and we take Ied = 119 + 8/n^. The ratio of the pressures 
behind the reverse and forward shocks is given by <1>ed = 0.39 — 0.6 exp (— n/4). The values 
for these parameters (both for s = 2) are fitted to the tabulations given in Chevalier (1982). 
The blast wave velocity is given by 



and the reverse shock velocity by 

Rr dRr I" 3 — s\ Vi, 

^ t dt — 3/ Ied 


(A6) 


(A7) 


The blast wave radius during the phase when the reverse shock is propagating through 
the ejecta core is given by (from equation 45 of Truelove & McKee 1999) 



(3-s) 



w 


3/2 



n 

3 (n - 3) 


A/2 

Jn 

n — 3 




(A8) 

where Wcore = f^corelf^ej and Vcore = y^l0/3-^(n — 5) / (n — 3) is the ejecta velocity at the 
core-envelope boundary and Vgj is the ejecta velocity at the outer edge of the envelope. In 
deriving equation A5 we have put Wcore — > 1- The ejecta structure function / (w) = fo for 
0 < 10 < Wcore and / (w) = /„/io” for Wcore ^ w < 1. Prom continuity at the core-envelope 
boundary fo = fn/w'^^ 3/47t as Wcore — >• 1- The time tcore where the reverse shock hits 

the ejecta core is determined from equation A5 by setting = lED’^caretcorei '^core 1 so 
that Vcore = Uej and solving for tcore with the result 


t - 

f 1 3(3-s)^ ' 

1 l/(3-s) 

1 

CO 

^core — ; — > 

^ED'^core | 


( 

lO 

1 

o 
1 


(A9) 


This generalizes equation 79 of Truelove & McKee (1999). Equation A8 is quite simple to 
work with for s = 0, but for s = 2 it requires the solution of a cubic equation in R-l^^, and 
for arbitrary s can be even more complicated. We adopt a simpler procedure of extending 
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the blast wave envelope solution into the core phase and matching it to the offset power law 
solution that is appropriate in the Sedov-Taylor limit, 


_n-3Rb _ d («-3)/2 

‘^b — j. ~ r. ^b 

n — s t o — s 


(AlO) 


where C = (5 — s) (10 — 3s) /Stt. The time at which these two solutions connect is derived 
by eliminating uj, between equations A6 and A9, and then substituting the envelope solution 
for equation A5 into the resulting expression for t and solving for t with solution 



2(n-s)/(5-n)(3-s) 


(All) 


The complete expression for the forward shock radius is 


Rb=< 


f (3 s) 


-I (5-s)/2(n-s) 


{tconn'^CMre) 


n— 3 


I (t>ED^^ri (n - 3) 

with the velocity given by equation AlO. 


2/(5-s) 




Stt 


(A12) 


The reverse shock trajectory through the ejecta core and into the Sedov-Taylor phase 
cannot in general be specified without recourse to numerical calculations. Truelove &; McKee 
(1999) show that for s = 0 SNRs the reverse shock velocity is approximately proportional 
to time and so introduce a constant acceleration parameterization, with the value of the 
acceleration being determined numerically. In general the acceleration turns out to be quite 
small, and the reverse shock velocity is almost constant. By the same methods it can be 
shown that for s = 2 SNRs the reverse shock velocity varies as \/f, and following reverse 
shock propagation into the ejecta core, we expect it to accelerate much less than in the s = 0 
case, and so during the core-Sedov-Taylor phase we hold Vr constant at the value given by 
equation A7 upon entry into the ejecta core. The reverse shock radius in this phase is given 


by 


Rr 


-I 


Rb — ^core) 
^ED^'core 


Vr In (^/^core) ^ 


(A13) 


To summarize, the blast wave radius and velocity are given by equations A45 and A5 for 
t < tconn, given by equation All, and are given by equations AlO and A12 for t > tconn- The 
reverse shock radius and velocity (with respect to the otherwise freely expanding ejecta) are 
given by equations A5 (divided by Ied) and A7 for t < tcore, with tcore given by equation A9. 
For t > tcore the reverse shock velocity is held constant and the radius is given by equation 
A13. 
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Truelove k, McKee (1999) give no guidance on the motion of the contact discontinuity. 
We use the results of Chevalier (1982) for the envelope phase, and during core propagation 
we assume that the contact discontinuity expands with 0.75 of the forward shock velocity. 
Detailed calculations indicate that the forward shock-contact discontinuity separation in- 
creases once the reverse shock has reached its maximum radius (Wang & Chevalier 2001). 
We estimate that this will occur once Cas A is about 770 years old for the n = 6 ejecta 
profile, the with projected time increasing rapidly with increasing n. An inconsistency arises 
with this assumption for the steeper models, in that the reverse shock radius can be ahead 
of the contact discontinuity. We suspect that our estimate of the reverse shock velocity is 
becoming inaccurate, and increase Vr in these cases so that Rr as given by equation A13 
is always less than Rc = 0.75i?&. This modification produces essentially no change in the 
spectroscopic parameters, Te and Uet in which we are interested. We have also tried several 
different analytic representations of the hydrodynamics and find our curves of Tg against Ugi 
to be insensitive to these differences. 


B. BLASPHEMER Simulations 


BLASPHEMER (BLASt Propagation in Highly EMitting EnviRonment) follows the 
time dependent ionization balance and temperatures of a Lagrangian plasma parcel as it 
passes through either the forward or reverse shock and then expands with the rest of the 
supernova remnant. In this appendix we repeat and update some of the description in 
Laming (2001b). Behind the reverse shock, the density Ug of ions with charge q is given by 


dvi 

[(^rr,^H-l “I” C^dr,g+l) (C^rr,g “1“ (^1) 


where Cion,q, Crr,q, Cdr,q ^re the rates for electron impact ionization, radiative recombination 
and dielectronic recombination respectively, out of the charge state q. These rates are the 
same as those used in the recent ionization balance calculations of Mazzotta et al. (1998), 
using subroutines kindly supplied by Dr P. Mazzotta (private communication 2000). The 
electron density Ug is determined from the condition that the plasma be electrically neutral. 
The ion and electron temperatures, Tj and Tg are given by 


dTi (Ti-Tg)E,9S/(9 + l) 

■5T = 


(B2) 


f®«=0 13 (ri-T,)E,A,/(?+l) T,dn. 2 dQ 

dt ' ATe^'^ J2q '^e dt SrigAlB dt ■ 


and 


(B3) 
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Here A is the atomic mass of the ions in the plasma. The last term dQ/dT represents 
plasma energy losses due to ionization and radiation. Radiation losses are taken from Sum- 
mers &; McWhirter (1979). At each time step Tj and Tg are modified by a further factor 
exp (— 4ueiAt/3r) and Ug and the n, by exp {— 2uexA</r). Here Vg^ is the expansion veloc- 
ity, which is assumed constant for all ejecta, and is given by the expansion velocity of the 
contact discontinuity discussed above. The radial compression/decompression of shells of 
shocked ejecta is treated by further modifying densities and temperatures by- Vgx/v'^j. and 
(vgx/v'gj.)^^^ respectively at each time step, where is Vgx at the previous time step. The 
plasma pressure is evolved according to adiabatic expansion in the same way, and the den- 
sities and temperatures further adjusted by {P/P'Y where a = 0.6 for densities and 0.4 for 
temperatures. This accounts for the compression of plasma due to radiation and ionization 
losses. The initial plasma pressure is P and the plasma pressure after losses is P', and the 
Lagrangian plasma element is recompressed adiabatically at each time step to restore its 
pressure to the adiabatic value after losses. Except for ejecta near the contact discontinuity, 
radiation losses are generally unimportant. 


Spitzer (1978) gives the timescale for an electron distribution to relax to a Maxwellian 


as 


^eq (®) c) — 


ZrnJ^ 


(ksTg) 


3/2 


(B4) 


In A 

where A is the so-called plasma parameter, the ratio of largest to smallest impact parameters 
for collisions. In supernova remnants In A ~ 40. The equilibration time for ions tgg (i, i) = 
tgq (e, e) y/mi/mg/Zf, and that for electron-ion equilibration is tgq (e, i) — tgq (e, e) rrii/mg/Zf 
where Zi is the ion charge. Accordingly we write 


dAT 

dt 


-o.iaz^rie 


AT 


(B5) 


which is equation (2), with AT = Ti — Tg. We consider a fully ionized gas with Ug = Zrii 
and 

^ m \ dTi dTe ^ 

— {niTi-\-ngTg) = ni— + ng— = 0. (B6) 

Solving these equations yields 


dt 

dt 


= 0.13 


Z^rig Tj-Tg 

z + 1 at |/^ 


-0.13 


Z^Tig Ti - Tg 

Z + 1 at!^^ ■ 


(B7) 

(B8) 


In deriving equations (14) and (15) these expressions are averaged over the ion charge* states 
in the plasma, and the expression for dTg/dt is modified by the inclusion of terms accounting 
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for the change in electron densitj<(l'iltetfcekJirBdMeii4.a^$’^aa?):^ici3«ij/(f4i^ and radiative and 
ionization losses, — (2/3neA:B) dQ/dt. Recombinations, which reduce the 
not result in an increase in the electron temperature in low density plasmas, since the energy 
of the recombined electron is radiated away, rather than being shared with the other plasma 
electrons as would be the case for three-body recombination in dense plasmas. 
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Fig. 1. — Locations of series of “Si” rich knots on northern limb. East is to the left and 
north is up. 





Fig. 2. — Detail of Cas A showing locations of “Si” rich knots at the jet base in the NE 
region. East is to the left and north is up. 
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Fig. 5. — Plots of Te against Ugi for varying ejecta envelope power laws. The solid lines 
give models with M^j = 2M©, E^i = 3 and prl = 21, for n = 5, 5, 6, 7, 9, and 12 from the 
top. The point at highest ngt for n = 5.5 and 6 corresponds to ejecta at the core-envelope 
boundary. For higher values of n this plasma undergoes thermal instability. The dotted lines 
indicate the variation of Tg with nj, for the modified models in Table 4 forn = 5.5, 6. Higher 
values of n change insignificantly in the modified models. Only fits with xVd-o.f. < 2 are 
plotted. 
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ionization age (cm"'*s) 


Fig. 6. — Plots of Te against ngt for varying ejecta envelope power laws. The solid lines give 
models with M^j = 2Mq, E 51 = 2 and pr^ = 14, for n = 5.5, 6, 7, 9, and 12 from the top. 
The point at highest ngt for n = 5.5, 6, and 7 corresponds to ejecta at the core-envelope 
boundary. For higher values of n this plasma undergoes thermal instability. The dotted lines 
indicate the variation of Tg with ngt for the modified models in Table 5 forn = 5.5, 6. Higher 
values of n change insignificantly in the modified models. Only fits with x^/d-O-f. < 2 are 
plotted. 
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Table 1. Fits to 0-Si rich Knots 


knot^ 

fcBTe (keV)^ 

net (cm“V) 

Si/O^ 

Pe/O 

riH (10^^ cm“^) 


X^/dof 

no. dof 

no. cts 

El 

1.78 

6.9e-fl0 

0.08 

0.05 

1.60 

233.1 

1.46 

160 

9832 


1.70-1.85 

6.0e+10-7.3e+10 

0.07-0.09 

0.04-0.06 

1.49-1.67 





E2 

1.40 

2.5e-hl0 

0.35 

0.04 

1.99 

199.9 

1.54 

130 

7976 


1.22-1.47 

2.2e-Hl0-2.7e-|-10 

0.31-0.37 

0.03-0.04 

1.96-2.05 





E3 

1.05 

3.2e-l-ll 

0.28 

0.17 

1.22 

236.7 

1.75 

135 

11897 


1.03-1.10 

2.5eH-ll-4.2e-fll 

0.24-0.39 

0.13-0.20 

1.16-1.28 





E4 

1.13 

1.8e+ll 

0.29 

0.07 

1.11 

184.7 

1.71 

108 

5994 


1.04-1.22 

1.5eH-ll-2.3e-|-ll 

0.27-0.35 

0.06-0.10 

1.03-1.24 





E5 

1.34 

2.3e+ll 

0,29 

0.12 

0.99 

236.9 

1.74 

136 

10797 


1.21-1.37 

2.2eH-ll-2.5e-Ml 

0.27-0.31 

0.09-0.12 

0.96-1.02 





E6 

1,23 

8.5e+10 

0.29 

0.05 

1.51 

271.8 

1.84 

148 

10452 


1.18-1.25 

8.2e+10-8.9e4-10 

0.28-0.34 

0.04-0.06 

1.49-1.58 





E7 

0.91 

2.0e+ll 

0.28 

0.04 

1.2 

332.8 

2.47 

135 

9389 

E8 

0.72 

4.2e-hll 

0.45 

0.04 

1.2 

296.0 

2.29 

129 

10316 

E9 

1.22 

3.1e+ll 

0.28 

0.13 

1.20 

163.6 

1.56 

105 

4924 


1.20-1.27 

1.8e-|-ll-3.6e-Hll 

0.23-0.36 

0.11-0.21 

1.17-1.28 





NNWl 

1.42 

l.le+11 

0.29 

0,12 

1.3 

335.8 

2.05 

164 

15126 

NNW2 

1.85 

1.2e-!-ll 

0.24 

0.11 

1.1 

760.7 

3.21 

237 

32844 

NNW3 

1,75 

1.4e-hll 

0.21 

0.09 

1.0 

548.1 

3.11 

176 

16483 

NNW4 

1.86 

5.1e+10 

0.32 

0.11 

1.49 

203.3 

1.42 

143 

8725 


1.63-1,95 

4.8e-|-10-5.9e4 10 

0.29-0.34 

0.10-0.12 

1.45-1.54 





NNW5 

1.49 

7.6e-l-10 

0.24 

0.08 

1.22 

242.9 

1.48 

164 

12960 


1.45-1.67 

6.9e-f-10-8.0e-l-10 

0.22-0.25 

0.07-0.09 

1.20-1.28 





NNW6 

1.10 

3.3e+10 

0.20 

0.00 

1.7 

305.0 

2.18 

140 

10584 

NNW7 

1.54 

1.7eH-ll 

0.18 

0.04 

0.84 

210.9 

1.70 

124 

6826 


1.44-1.62 

1.6eH-ll-1.9e411 

0.17-0.19 

0.03-0.04 

0.76-0.86 





NEl 

1.11 

1.6e-|-ll 

0.06 

0.08 

1.08 

87.4 

1.00 

87 

4124 


0.97-1.19 

1.2e-|-ll-2.7e-|-ll 

0.05-0.08 

0.06-0.10 

0.95-1.24 





NE2 

1.12 

6.5eH-10 

0.40 

0.04 

1.24 

189.1 

1.70 

111 

6009 


1.03-1.22 

6.2e-|-10-6.8e410 

0.38-0.52 

0.03-0.04 

1.10-1.29 





NE3 

1.54 

3.2e-|-10 

0.29 

0.08 

1.35 

119.6 

1.34 

90 

4155 


1.44-1.61 

3.0e+10-3.4e-hl0 

0.27-0.31 

0.07-0.09 

1.32-1.37 





NE4 

1.53 

3.0e+10 

0.60 

0.04 

1.29 

266.8 

1.89 

141 

10884 


1.46-1.59 

2.8e-|-10-3.1e-|-10 

0.57-0.62 

0.03-0.04 

1.27-1.33 





NE5 

1.51 

3.3e-l-10 

0.44 

0.03 

1.4 

284.1 

2.47 

115 

6961 

NE6 

1.32 

2.8e+10 

0.34 

0.00 

1.61 

191.4 

1.48 

129 

7578 


1.29-1.46 

2.3e+10-2.9e-l-10 

0.32-0.36 

0.00-0.003 

1.60-1.63 





NE7 

1.67 

4.5e-f-10 

0.25 

0.06 

1.54 

280.4 

1.76 

159 

10016 


1.62-1.73 

4.3e+10-4.7e+10 

0.24-0.26 

0.05-0.06 

1.52-1.56 





NE8 

1.90 

4.3e-HlO 

0.26 

0.04 

1.52 

248.0 

1.53 

162 

9776 


1.67-1.98 

4.3e-l-10-4.8eH-10 

0.23-0.27 

0.03-0.04 

1.48-1.56 





NE9 

1.51 

6.1e+10 

0.21 

0.05 

1.31 

308.2 

1.78 

173 

14081 


1.47-1.59 

5.7e-hl0-6.4e+10 

0.19-0.22 

0.04-0.05 

1.28-1.32 





NEIO 

1.67 

9.1e-f-10 

0.06 

0.09 

1.31 

221.9 

1.66 

134 

8451 


1.53-1.75 

8.2e-K0-l.le+ll 

0.05-0.08 

0.07-0.10 

1.27-1.38 





NEll 

1.10 

9.1e+10 

0.12 

0.17 

1.32 

162.7 

1.75 

93 

4812 
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Table 1 — Continued 


knot^ (keV)^ 

7let (cm ^S) 

Si/O' 

Fe/O 

TiH (10^^ cm~^) 

X^/dof no. dof no. cts 

0.90-1.26 

6.7e-|-10-1.3e-fll 

0.09-0.15 

0.14-0.21 

1.22-1.41 



^Knots are numbered in the E and NNW regions starting with the innermost and working out to the limb. The knot numbering 
in the NNE region is given in Figure 2. 

^Electron temperature in keV. 

*^Element abundance ratio by number relative to solar values of Anders & Grevesse (1989). Note that these are superseded 
by Grevesse & Sauval (1998), Allende Prieto, Lambert, h Asplund (2001), and Allende Prieto, Lambert, & Asplund (2002). In 
particular abundances relative to O increase by 1.75. 


Table 2. Cas A Ejecta Profile Models Mej = 2Mq, E^\ = 3, fyrl = 21 


n 

Vf, (320yrs) 
km s"^ 

rb (320yrs) 
pc 


Tr (320yrs) 
pc 

hItt 

Vcore^ 
km S~^ 

t ^ 

•'Core 

yrs 

t ^ 

iconn 

yrs 

^rad^ 

yrs 

^rad^ 

Mo 

5.5 

4044 

1.85 

0.72 

0.92 

2.00 

7098 

39.1 

50586 

_ 

0 

6 

4942 

2.15 

0.75 

0.97 

2.21 

9163 

21.5 

632 

- 

0 

7 

5139 

2.36 

0.71 

1.13 

2.08 

11223 

10.1 

120.2 

10.05-10.25 

0.012 

8 

5108 

2.38 

0.70 

1.35 

1.77 

12294 

5.99 

72.4 

4.55-10.4 

0.41 

9 

5095 

2.40 

0.70 

1.58 

1.52 

12959 

3.97 

55.0 

2.67-11.1 

0.61 

10 

5088 

2.41 

0.69 

1.79 

1.34 

13414 

2.84 

45.7 

1.76-11.4 

0.69 

11 

5083 

2.41 

0.69 

1.99 

1.21 

13745 

2.14 

39.8 

1.25-11.9 

0.75 

12 

5079 

2.41 

0.69 

2.17 

1.11 

13997 

1.67 

35.6 

0.93-12.5 

0.79 


®^Forward shock expansion parameter. 

^Free expansion velocity of ejecta core-envelope boundary. 

^Time following explosion when reverse shock enters ejecta core. 

^Time when blast wave solutions are connected. 

®Time interval for which ejecta passing through the reverse shock cools to optically emitting temperatures within 320 
years. 

^Mass of gas that can cool to optically emitting temperatures within 320 years of explosion. 


Table 3. Cas A Ejecta Profile Models Mej = 2Mq, E^i = 2, prl = 14 


n 

VI (320yrs) 
km s“^ 

Tb (320yrs) 
pc 

»?* 

Tr (320yrs) 
pc 

n/fr 

Vcore^ 
km s“^ 

# c 

icore 

yrs 

/ d 

t'conn 

yrs 

trad* 

yrs 

^rad^ 

Mq 

5.5 

3928 

1.79 

0.72 

1.08 

1.66 

5795 

71.8 

92930 

_ 

0 

6 

4698 

2.04 

0.75 

1.17 

1.75 

7482 

39.5 

1162 

- 

0 

7 

5239 

2.27 

0.76 

1.29 

1.76 

9163 

18.6 

221 

- 

0 

8 

5177 

2.32 

0.73 

1.44 

1.62 

10038 

11.0 

133 

9.25-15.5 

0.26 

9 

5153 

2.35 

0.72 

1.59 

1.48 

10581 

7.3 

101 

5.35-16.5 

0.50 

10 

5139 

2.36 

0.71 

1.74 

1.36 

10953 

5.2 

84.0 

3.45-17 

0.60 

11 

5129 

2.37 

0.71 

1.87 

1.27 

11223 

3.93 

73.1 

2.45-17.5 

0.66 

12 

5121 

2.37 

0.71 

2.00 

1.19 

11429 

3.07 

65.5 

1.815-18 

0.70 


^'Forward shock expansion parameter. 

^Free expansion velocity of ejecta core-envelope boundary. 

*^Time following explosion when reverse shock enters ejecta core. 

^Time when blast wave solutions are connected. 

®Time interval for which ejecta passing through the reverse shock cools to optically emitting temperatures within 
320 years. 

^Mass of gas that can cool to optically emitting temperatures within 320 years of explosion. 


Table 4. Models = 2.40 pc, pcore = l-22e6 g cm ^s^, and prl = 21 H atom cm ^ pc 


n 

Mej 

Mq 

£^51 

10^^ ergs 

Vb (320yrs) 
km s“^ 

Tb (320yrs) 
pc 


Tr (320yrs) 
pc 

n/rr 

Mtj (n - 3) /n/ufore** 
g cm”^s® 

5.5 

1.72 

5.0 

5236 

2.391 

0.716 

0.832 

2.87 

1.226e6 

6 

1.72 

3.85 

5530 

2.405 

0.752 

0.867 

2.77 

1.219e6 

7 

1.85 

3.10 

5181 

2.396 

0.707 

1.09 

2.20 

1.226e6 

8 

1.95 

3.05 

5132 

2.403 

0.699 

1.35 

1.78 

1.225e6 

9 

2 

3 

5095 

2.398 

0.695 

1.58 

1.52 

1.219e6 

10 

2.04 

2.99 

5084 

2.401 

0.693 

1.79 

1,34 

1.218e6 

11 

2.07 

2.98 

5075 

2.402 

0.691 

1.96 

1.23 

1.226e6 

12 

2.08 

2.97 

5065 

2.402 

0.690 

2,14 

1.12 

1.218e6 


^Forward shock expansion parameter. 
*^Ejecta core density in velocity space. 



Table 5. Models rt = 2.35 pc, pcore = 2.25e6 g cm and = 14 H atom cm ^ pc 


n 

Mej 

Mo 

£^51 

10®^ ergs 

Vb (320yrs) 
km s“^ 

Tb (320yrs) 
pc 


Tr (320yrs) 
pc 

n/rr 

M^j (n - 3) /n/v^ore'' 
g cm“^s^ 

5.5 

1.815 

4.15 

5134 

2.345 

0.716 

1.28 

1.92 

2.146e6 

6 

1.815 

2.8 

5395 

2.347 

0.752 

1.21 

1.94 

2.248e6 

7 

1.875 

2.15 

5338 

2.349 

0.744 

1.29 

1.82 

2.256e6 

8 

1.95 

2.05 

5212 

2.349 

0.726 

1.44 

1.63 

2.223e6 

9 

2 

2 

5153 

2.345 

0.719 

1.59 

1.48 

2.240e6 

10 

2.035 

1.985 

5130 

2.348 

0.715 

1.73 

1.36 

2.238e6 

11 

2.065 

1.965 

5105 

2.347 

0.712 

1.85 

1.27 

2.267e6 

12 

2.075 

1.95 

5085 

2.347 

0.709 

1.96 

1.20 

2.267e6 


^-Forward shock expansion parameter. 
^Ejecta core density in velocity space. 


